Polymers grafted to porous membranes 
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We study a single flexible chain molecule grafted to a membrane which has pores of size slightly 
larger than the monomer size. On both sides of the membrane there is the same solvent. When 
this solvent is good, i.e. when the polymer is described by a self avoiding walk, it can fairly easily 
penetrate the membrane, so that the average number of membrane crossings tends, for chain length 
N — > oo, to a positive constant. The average numbers of monomers on either side of the membrane 
diverges in this limit, although their ratio becomes infinite. For a poor solvent, in contrast, the 
entire polymer is located, for large N, on one side of the membrane. For good and for theta solvents 
(ideal polymers) we find scaling laws, whose exponents can in the latter case be easily understood 
from the behaviour of random walks. 



Long chain polymers in dilute solutions have impor- 
tant technical applications, but they also form an in- 
teresting model for studying scaling and renormalization 
group techniques Q, 0,13 ■ As a consequence, most of 
their properties have been studied in detail and are well 
understood, at least as far as equilibrium is concerned. 
In this paper we want to study a new effect which, to 
our knowledge, has not been considered in the previous 
literature, although also it has a number of possible ap- 
plications. 

Physically, we ask ourself what is the influence of the 
solvent on the ability of a polymer to penetrate a porous 
membrane. If the pores are very big, typically of sizes 
larger than a collapsed globule but smaller than the Flory 
diameter of a random coil in a good solvent, the collapse 
might help in penetrating them. On the other hand, 
we expect the opposite effect for small pores, typically 
just slightly larger than a single monomer. In that case a 
polymer in good solvent can go through a pore by having 
one monomer pass after the other |4|. For a collapsed 
polymer this should be less easy, since monomers tend 
to cluster together. But, as far as we know, no detailed 
studies of this exist. 

Technically, we consider self avoiding walks (SAWs) of 
N steps on a simple cubic lattice. In a good solvent at 
temperature T — > oo, the only monomer-monomer inter- 
action is the excluded volume interaction. In a less good 
solvent, we include an attractive monomer-monomer con- 
tact potential which acts between unbonded monomer 
pairs on nearest neighbour sites on the lattice, and has 
strength — e, giving for each such contact a Boltzmann 
factor q = e e l KbT . The theta point for this model is at 
qe = 1.3087 ± 0.0003 [|. The membrane is located in 
the plane z = 0. Its pores are of unit size, and are ar- 
ranged in checkerboard fashion: sites with even x + y and 
z = arc forbidden for the walk, while all other sites are 
allowed. The walks start at the origin x = y = z = 0, 
i.e. one end is grafted to one of the forbidden sites of the 
membrane. 



We simulate this model using the pruned-enriched 
Rosenbluth method (PERM) |j. This is the most ef- 
ficient known method for simulating theta polymers, but 
it is also very efficient for conditions not too far from the 
theta point, including athermal SAWs. For the latter, 
the pivot algorithm is even more efficient if the SAWs 
are not constraint by obstacles like our membrane. If 
there are many such obstacles, PERM becomes more ef- 
ficient. We have not checked whether the pivot method 
would be more efficient for the present problem, but this 
is of little relevance since we were anyhow able to obtain 
very high statistics data for very long chains (N — 16000 
at the theta point, N = 10000 for athermal SAWs). 

To measure how well the polymer can penetrate 
through the membrane, we use two different "order pa- 
rameters" (i.e., observables). The first is simply the num- 
ber Nq of points on the walk in the plane z — (excluding 
the starting point). For the other, we first define N + as 
the number of monomers in the upper half space z > 0, 
and iV_ as the number of monomers at z < 0. Obviously, 
N + +N- + No = N. The second order parameter is then 

a= NfTW (1) 

In Fig. 1 we plot a against N, for six different values 
of q. Apart from the strong finite size effects at small 
N (a vanishes trivially for N < 4), there seems to be a 
power behaviour for q < q$. Indeed, when plotting a^/N 
against N (Fig. 2), we see that 

a ~ N- 1/2 (2) 

for theta polymers, while a seems to decrease with a 
smaller power for polymers in a good solvent. Due to 
the aforementioned scaling corrections it is hard to give 
a more precise estimate, our best fit is a ~ jy-o.38 f or 
SAWs. For collapsed chains a decreases much faster, but 
it is not clear whether it will finally follow a power law 
(with power > 1) or whether it falls faster. 
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FIG. 1: Log-log plot of a against N, for six different values 
of q. Data decrease with q. Statistical errors are not larger 
than the thickness of the lines. 
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FIG. 2: Same data as in Fig. 1, but multiplied by N 1/2 . 
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FIG. 3: Log-linear plot of (iVo) against TV, for the same six 
different values of q. Again, data decrease with q. The rapid 
decrease of the data for q = 1.37 and q = 1.44 at very large 
TV is presumably due to statistical fluctuations. 



The average No is plotted against N, for the same six 
values of q, in Fig. 3. We see huge even/odd oscillations 
for small N. Apart from that, the results are somewhat 
complementary to those in Figs. 1 and 2. For q = qg we 
see that Nq converges quickly to a constant, 



0.5703 ± 0.0003 



for N 



oo 



(3) 



Again scaling corrections are more important for good 
solvent conditions. But although Nq keeps increasing 
with N in this case, it seems to converge nevertheless to 
a constant. This is definitely no longer true for q > qg. 
There, Nq seems to converge to zero. We cannot exclude 
that this happens according to power laws, N Q ~ N~P, 
but then the exponent would be non-universal, and 
would increase strongly with q. We consider it much 
more likely that AT decreases, in the collapsed phase, 
faster than any power of N. 

Let us now try to understand these results theoreti- 
cally, and let us start with q = qg. For simplicity we will 
neglect the logarithmic corrections holding at the theta 
point d H 0, H, anc ^ describe the polymers by ideal 
random walks. In this case, the walks are Markovian. 
Whenever such a walk arrives at z = 0, it will go in the 
next step to z > and z < with equal probability, 
and both possibilities will evolve later in the same way 
(except for mirror symmetry). Thus, as far as scaling is 
concerned, this is the same situation as for an impenetra- 
ble wall. In this case |9( the monomer density decreases 
near the wall as z 2 , and the number of monomers just 
touching the wall goes to a constant for N — > oo. This is 
in agreement with Fig. 3. 

The probability that the n-th monomer, with -C n <C 
N, is at distance z from the wall scales more precisely as 



7n 



3/2 



where the n-dependence is forced by the fact 



that (z) ~ \pn. Thus the probability that the wall is 
touched by the n-th monomer scales as n~ 3 / 2 , and this 
will also be the scaling of the probability that the walk 
passes through the membrane at "time" n. From this 
and the Markov property follows immediately Eq.JSJ. 

These simple arguments break down when q ^ qg. 
Now, it does matter that a walk which passes through a 
pore ends up on the opposite side of the membrane. For 
q < qg (good solvent conditions) it feels then a weaker 
excluded volume effect, as compared to the case where it 
stayed on the same side of the membrane. Thus passing 
through the membrane is (entropically) favoured, leading 
to larger values of Nq and a. On the other hand, pass- 
ing through the membrane is (energetically) disfavoured 
for bad solvent conditions, and Nq and a are reduced in 
this case. Apart from this very qualitative argument we 
have not been able to give a more thorough theoretical 
discussion. For SAWs grafted to a hard wall, the average 
number of monomers touching the wall goes to a con- 
stant for N — > oo 0, [n| . Thus it is not too surprising 
that the same holds also when the wall is replaced by a 
porous membrane, and only monomers within pores are 
counted. But it is not too obvious either, as the same 
analogy does not hold for collapsed polymers. While we 
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FIG. 4: Number of monomers at distance z — ±1 from the 
membrane, Ni, plotted against IniV for several values of q. 
While Ni converges to a constant when q < q$, it seems to 
diverge with N for q > qg. 
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FIG. 5: Logarithm of the partition sum divided by ix and 
multiplied by N 1 11 , plotted against In N. The values used 
for fi and 71 are 4.684038 and 0.679. They were fixed by de- 
manding that the leading corrections to scaling go like iV -1 / 2 . 



found here that N goes to zero for TV — > 00, the number 
Ni of monomers touching the wall or membrane (z = ±1; 
here it does not matter whether the membrane has pores 
or not) diverges with N, when the solvent is bad (q > qg; 
see Fig. 4). 

Finally, let us discuss how the porous membrane effects 
the partition function and the density of states g{E), 
where E = —em with m being the number of monomer- 
monomer contacts. For SAWs grafted to an impenetrable 
wall, the partition sum scales as |lfj| 



where /i = 4.684043(10) is the same effective connectivity 
as for ordinary SAWs, but 71 is a modified critical expo- 
nent, 71 = 0.679(2). Numerically (see Fig. 5) we found 
-2.8 
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FIG. 6: Ratio between the density of states for polymers 
grafted to a porous membrane, divided by the analogous den- 
sity of states for free polymers without a membrane. The 
logarithm of this ratio is plotted versus —E/N, for two dif- 
ferent values of N. 



that the same ansatz describes also the large TV behaviour 
of athermal SAWs grafted to the porous membrane (in 
Fig. 5 we used our present best estimates 71 = 0.679(1) 
and ix — 4.684038(6)). This is not trivial, since we had 
just argued that the possibility of penetrating the mem- 
brane should increase the partition sum. It does so in- 
deed, but this does not affect the asymptotic behaviour. 

Instead of showing the density of states g(E) itself, de- 
fined as the number of configurations with energy —E, we 
show in Fig. 6 the ratio g(E)/go(E), where go(E) refers 
to polymers in absence of a membrane. More precisely, 
we plot the logarithm of this ratio against —E/N for two 
values of N. Here we have assumed that e = 1, i.e. each 
contact contributes exactly one unit of energy. We see 
that g(E)/go(E) decreases with the absolute value of E, 
i.e. the effect of the membrane increases with \E\. This 
is of course in agreement with our other measurements. 
A more careful inspection shows that the curves in Fig. 6 
are not straight, but are cup-concave around E — 0.47V. 
This means that g{E) is less convex than go(E) in the 
energy range which gives the dominant contributions to 
the partition sum near the theta collapse, i.e. the theta 
transition becomes "harder" by the presence of the mem- 
brane. For the present model this is a rather weak effect, 
but is is not inconceivable that it is much stronger in some 
other microscopic realization. In such a model, grafting 
polymers to porous membranes may make the collapse 
transition look first order for intermediate values of N. 

We thank Walter Nadler for very helpful discussions. 
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